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I. INTRODUCTION 



This thesis is but a part of an ongoing research project to develop new sonar 
signal processing algorithms capable of quickly and accurately solving target lo- 
calization problems. Present technology and doctrine dictate that several lines of 
bearing to a target be obtained before a sonar fire control solution can be com- 
puted. Obtaining these lines of bearing is a time consuming and often dangerous 
task due to the increased probability of counterdetection and, as a result, evasive 
maneuvering and defensive action on the part of the target. A sonar system ca- 
pable of providing timely, accurate target localization while minimizing own ship 
maneuvering would result in longer firing ranges and, therefore, a reduction of the 
threat to one’s own ship. 

Several recent papers [Refs. 1-3] have discussed the application of a complex 
least-mean-square (LMS) adaptive algorithm [Ref. 4] to bearing estimation prob- 
lems using a linear array of sensors. When a linear array is used, only a bearing 
angle to the source can be estimated, an estimate of the depression angle 6 ^ can- 
not (Figure 1.1). While an estimate of the bearing angle is useful in the localization 
problem, a better tool would be an algorithm which provides estimates of both the 
bearing angle and the depression angle. A frequency domain adaptive beamform- 
ing algorithm for planar arrays that solves multiple, broadband target localization 
problems and provides estimates of both the bearing angles and the depression 
angles from the center of the planar array to the targets has been developed and 
initially tested by Ziomek and Chan [Ref. 5j. 
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Figure 1.1. Target Location Geometry. 



2 



The research performed in this thesis was to continue the work of Ziomek 
and Chan [Refs. 5, 6] and fully evaluate the capabilities of the complex, least- 
mean-square, frequency domain, adaptive beamforming algorithm they developed 
via computer simulation studies. The computer simulation studies were designed 
to test the algorithm’s multiple broadband target localization capability, its full 
angular coverage capability, and its angular resolution as a function of the input 
signal- to- noise ratio (SNR) at a single element in the array, sampling rate, harmonic 
number, and the number of iterations of the algorithm. 

Each target was modeled as a broadband sound source. As a result, the 
frequency spectrum of the output signal from each element of the planar array 
contains several frequency components. An estimate of the bearing and depression 
angle for each frequency component is provided as a result of processing the output 
frequency spectrum from each element in the array through the frequency domain 
adaptive beamforming algorithm. Therefore, if each target exhibits at least one 
unique frequency component (or spectral line), then all targets can be located. 

Full angulair coverage is the ability to localize a target regardless of its relative 
position to the array. The broadside case is the easiest since it is at this position 
that the far-field beam pattern beamwidth is its narrowest. The endfire case is the 
most difficult since the far-field beam pattern beamwidth is the broadest at this 
point. The full angular coverage and multiple broadband target capabilities were 
tested simultaneously. 

Baseline results were the first assembled. Baseline results are defined as beeir- 
ing and depression angle estimation errors (measured in degrees) as a function of 
sampling rate, harmonic number, and the number of iterations of the algorithm 
for the “no noise” case. Identical cases were run using additive, zero mean, white. 
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guassian noise to corrupt the output signals from each element of the planar ar- 
ray. Average bearing and depression angle estimation errors were then plotted as 
a function of the input SNR at a single element of the planar array, the sampling 
rate, and the harmonic number. 

Chapter II describes the theory used in the development of the frequency do- 
main adaptive beamforming algorithm. The construction of the direction cosine 
estimates, the angle estimates, and phase “unwrapping”, integral parts of the al- 
gorithm, will be presented in detail. 

Chapter III contains computer simulation results and an explanation of these 
results. The results which are presented consist of graphical representations of 
the average estimation errors of the bearing and depression angles at two distinct 
signal-to-noise ratio levels for four cases. These cases include targets evaluated 
at broadside, endfire, random placement, and tairgets which share a spectral line. 
Conclusions concerning the effect of harmonic number, sampling rate, number of 
iterations of the algorithm, and SNR are made. The Appendix contains tabular 
numerical data for all results. 

Chapter IV will present final conclusions and recommendations for further 
research. 

Results of this thesis research work have been incorporated into a research 
paper by Ziomek and Behrle which was presented at the Twenty-First Annual 
Asilomax Conference on Signals, Systems, and Computers [Ref. 7]. 



4 



II. THEORY 



This chapter is designed to present a brief, yet comprehensive, development of 
the frequency domain adaptive beamforming algorithm which is evaluated in this 
thesis. The goal of the algorithm is to process the output electrical signals from 
a planar array of sensors and provide estimates of direction (both bearing and 
depression angles) and frequency content of the acoustic fields incident upon the 
array. The development of the algorithm presented here is similar to the analysis 
section of a paper by Ziomek and Behrle [Ref. 7], but this chapter provides a more 
detailed explanation for several parts of the algorithm. A more in-depth analysis 
of the algorithm is provided by Ziomek and Chan [Ref. 5] and Chan [Ref. 6]. 

The frequency domain adaptive beamforming algorithm is based on using the 
complex frequency domain data R{q,m,n) from all M X AT elements of a planax 
array. The objective of the adaptive filter used in the algorithm is for the filter 
to converge to a set of phase weights such that the array output signal will match 
a reference signal. From this set of phase weights, estimates of the bearing and 
depression angles as a function of the harmonic number q can be made. 

Consider a. M x N plcinar array of point source elements lying in the XY 
plane (Figure 2.1). These elements are equally spaced where d* ^md dy are the 
interelement spacings in the X and Y directions, respectively, and M and N are 
the total number (odd) of elements in the X and Y directions, respectively. The 
random bziseband output electrical signal, r(£, m, n), is composed of a deterministic 
signal, y(£, m, n), and a random receiver noise component, n(£, m,n). The output 
signal at time instant £ and element (m, n) in the aurray is given by 
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Figure 2.1. Planar Array Geometry. 
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r(£, m,n) = y{t,m,n) + n(£, m,n), L = —L \ ... ,0, ... L' 



m = (2.1) 
n = 

where 

L' = {L- l)/2, 

M* = {M- l)/2, 

N' = {N -l)/2, 

and 

L > 2/C + 1 (2.5) 

is the total number of time samples that must be taken per element to avoid aliasing 
when the deterministic signal y{£,m,n) is composed of K harmonics [Ref. 8; p. 
164]. The minimum value for L (i.e., L = 2K + 1) is obtained by sampling at 
the Nyquist Rate. If To seconds is the length of the data record (or fundamental 
period) recorded at each element of the cirray, then the fundamental frequency (or 
the FFT bin spacing) is given by 



( 2 . 2 ) 

(2.3) 

(2.4) 



/o = l/To Hz. (2.6) 

Therefore, the highest frequency which can be contained within the deterministic 
part of the received signal is 

/max = Kfo Hz (2.7) 

and the minimum sampling frequency of the received output electrical signal r(£, m, n) 
is 

fs = LjTo Samples/sec (2.8) 
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where L must satisfy equation 2.5 [Ref. 8: p. 164]. 

To obtain the complex frequency domain samples of the received signal, the 
discrete Fourier transform (DFT) with respect to the time index I of equation (2.1) 
is taken. This action yields 

R{q, m,n) = Y (q, m, n) + N(q, m, n), q = -L', . . . , 0, . . . , L' 

m = (2.9) 
n = -iV',...,0,...,iV' 

where the index q represents the harmonic number. 

Consider a single general plane- wave field, g -b , propagating in the 

±no direction, incident upon the planar array as shown in Figure 2.2. If ^ is a 
unit vector, g{t) an arbitrary baseband function, and c the speed of sound in the 
medium measured in meters per second, then the deterministic part of the output 
electrical signal at time instant £ and element (m, n) in the array is given by [Ref. 
8: p. 160] 

rfi ^ f^rr - uomdx +vondY\ 
y{e,m,n) = g itTs + I (2.10) 

where Ts is the sampling period in seconds, and uq and uq are dimensionless 
direction cosines with respect to the X and Y axes, respectively. 

Upon taking the DFT of equation 2.10, we obtain the corresponding frequency 
spectrum given by [Ref. 8: pp. 162-166] 

Y{q,m,n) = Lcgexp{+ j2iTqfoUomdx/c)exp{+ j2TrqfoVondY /c) (2.11) 

where Cg, q = are the complex Fourier series coefficients that 

can be used to represent the baseband function g{t) by a finite Fourier series with 
K harmonics during the time interval jt| > Tq/ 2. The coefficients are given by 

c, = Y{q, 0, 0)/L, q = -L\ . . . , 0, . . . , i'. (2.12) 
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Figure 2.2. General Plane-Wave Field g Propagating in the 

direction. 
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The dimensionless direction cosines, Uq and vq, are given by 



uo = sin 9o cos V’o 
vq = sin ^0 sin t/>o 



(2.13) 

(2.14) 



where is the depression angle and is the bearing angle. These angles are not 
known a priori. 

If we consider multiple plane waves incident upon the planar array, equation 
(2.10) can be generalized and the deterministic part of the output electrical signal 
at time instant C and element (m,n) is 

^ f^rr. - uokmdx +vokndY\ 

y{e,m,n) = 2^gk yTs + ^ J (2.15) 

where gk{t), uok ', and vok are the arbitrary beiseband function and direction cosines 
associated with the kth sound source. 

A brief description of the development of the frequency domain adaptive beam- 
forming algorithm is presented next. The discussion that follows is based, in part, 
on the material in Ref. 7. 

Define the complex estimation error at harmonic q as 



4q) = - Kq)- 

The reference signal s{q) is defined as 

^ M' N' 

“ r.MM E E m, n)| exp [+jlR(q, 0, 0)] 



(2.16) 



and 



LMN 

^ M’ N’ 

~ T.MN E E c( 5 , m)d(q, n)R{q, m, n) 



(2.17) 



(2.18) 



m = -M* n=-iV' 



= c^{<t)R(<t)d(q)KLMN) 



10 



is the estimate of s{q), where c(q,m) and d{q,n) are the unit magnitude complex 
weights in the .Y and Y directions, respectively, c{q) and d(q) are the M x 1 and 
X 1 complex weight vectors in the X and Y directions, respectively, which are 
given by 

c(?) = [c{q,-M'), . . . ,c(q,0), . . . ,c(q, M')f (2.19) 

and 

d(g) = ld(q. -N')... ,,d(q, 0), . , . . d(g, N'f, (2.20) 



and R(q) is the M x N complex data matrix given by 

R{q,-M\-N') ... R{q,-M',0) ... R{q,-M\N'Y 



&(g) = 



m.g,0,-N') 



R{q,0,0) 



R(g,0,N') 



, R(q,M',-N') ... R(g,M',0) ... R(q,M',N') 

Now define the (M + N) x 1 complex weight vector w(q) as follows: 



w(q) = 



.d(g) 



( 2 . 21 ) 



( 2 . 22 ) 



We can obtain the original weight vector in the X direction via 



c(q) = A w(q) 



(2.23) 



where 



A = 



/ 1 


1 Q 


M xM 


MxN 



(2.24) 



is a M X (M + N) matrix, / is a M x M identity matrix, and 0 is a M x iV null 
matrix. Similarly, we can obtain the original weight vector in the Y direction via 



d(q) = B w(q) 

where 



(2.25) 



11 



(2.26) 



B = 



0 

NxM 



I 

NxN 



is a. N X [M + N) matrix, 0 is a iV x Af null matrix, zmd I_\s a N x N identity 
matrix as indicated. 

Substituting equations (2.23) and (2.25) into equation (2.18) yields 



5 ( 9 ) = m!'{q)Z{q)w{q)l{LMN) (2.27) 



where 

Z{q) = A^R{q)B (2.28) 

is a (A/+iV) X (M+N) complex matrix. The complex weight vector that minimizes 
the mean-square error £^{|e(g)p} is given by [Ref. 5, Ref. 7] 

Wi+i{q) = ULi{q) + 2/x,e.(<7) [Z{q) + Z^{q)]* Wi{q),i = 0,1,2,... (2.29) 

where 

ei{q) = s{q) - Si{q) (2.30) 

is the estimation error after the iteration, 

s»(9) = wJ{q)Z{q)wi{q)/{LMN) (2.31) 

is the estimate of s{q) after the iteration, and the step-size parameter m is 
given by 

M, = Mo = {<7l+ol)-\ t = 0,1,2,... (2.32) 

where /xq is a constant, and is equal to the inverse of the sum of the signal and 
noise power, <Xy and < 7 ^, respectively, at the center element of the array. After each 
iteration of the algorithm, each component of the complex weight vector ( 9 ) 
is normalized by its respective magnitude in order to maintain unit magnitude. 
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Once the complex weight vector, w-_^^{q), converges to a steady-state value, 
w^g{q), the steady-state complex weight vectors c^^iq) and d^^iq) can be obtained 
from equations (2.23) and (2.25), respectively. The estimates of the depression and 
bearing angles at harmonic q, 6o{q) and i^o{q), respectively, are given by [Ref. 8: 
p. 175] 



and 



^0(9) = sin 



-1 



+ 9#0 ( 2 . 33 ) 



^ 0 ( 9 ) = tan ‘ [Ci"*(9)/ao*(9)] , 9 # 0 



(2.34) 



where UQ^(q) and Vo^(q) represent estimates of the direction cosines obtained by 
using a least-squares fit to the “unwrapped” steady-state phase weights d'^^(q,m) 
and 0^^{q,n), respectively. In the absence of noise 

d'^,,{q,m) = ±2irqfoUo{q)mdx/c, = -M' '^0 

and 

= ±2TrqfoVo{q)ndY/c, ^ = -N> '^0 



It should be noted that if r(C,m,n), y(£,m,n), and n(i,m,n) of equation 
(2.1) are baseband complex envelopes, then qfo must be replaced by (/c -|- qfo) in 
equations (2.35) and (2.36), where fc is the carrier frequency in hertz. In addition, 
equations (2.33) and (2.34) must now be evaluated at 5 = 0. Substitution of 
equations (2.17) and (2.18) into equation (2.16) yields the following expression for 
the steady-state estimation error 

^ M’ N‘ 

m=^M> n=-N> 

-exp {+j{ei{q,m) + <i>Z{q,n) + lR{q,m,n)\] (2.37) 



exp[-|-;Zi?(9, 0,0)1 



13 



where 



R{q,m,n) = m, n)| exp[+; m, n)], 

Css{q,m) = a,,{q,m)exp[+je^^{q,m)] = exp [+jO'^^{q,m)\ , 



(2.38) 

(2.39) 



and 

dss{q,n) = b,s(q,n)exp[+j(f>'^^{q,n)] = exp [+;>^,( 9 , n)] (2.40) 

where a, ,( 5 , 771 ) = 1 and 65 ,( 5 , n) = 1 real, unit magnitude, amplitude weights 
and ^^( 5 , 771 ) and <^^,( 5 ,n) are real, “wrapped”, phase weights. 

To ensure obtaining the correct depression and bearing angle estimates, ^ 0 ( 9 ) 
and V’o(?)) respectively, the steady-state phase weights need to be “unwrapped.” 
The “unwrapped” steady-state phase weights can take on values outside the closed 
interval [— tt, tt] and therefore ensures full angular coverage capability (i.e., 0 < 
^ 0 ( 9 ) < 7t/ 2 and 0 < tpo{q) ^ Stt). The use of “wrapped” phase weights is necessi- 
tated by computer programming limitations. 

The question now arises how to obtain the “unwrapped” steady-state phase 
weights in order to obtain least-squares estimates of the direction cosines at each 
harmonic. The first step is to force the “wrapped” steady-state phase weights, 
^^,( 5 , 777 ) and <^^,( 9 ,n), to be equal to zero at the center of the planar array (i.e., 
777 = 0,71 = 0). This is accomplished by multiplying each component of c,,( 5 ) 
and d^giq) (equations (2.19) and (2.20)) by exp[— ^^^^',( 5 , 0)] and exp[— j(; 6 J“,( 5 , 0)], 
respectively. 

The next step is to obtain rough estimates of the direction cosines, uo{q) and 
Vo{q). The “wrapped” steady-state phase weight, 0'l^^(q,m), can be expressed as 
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i = 0, ±1, ±2, . . . 



±\2TtqfoUo{q)mdilc] + 2iT:, g = -L', . . . ,0, . . . ,L' (2.41) 

where i is chosen to ensure that the value of 0f^{q, m) is within the interval [— tt, tt). 
In the computer simulations presented in this thesis, the interelement spacing is 

di = dy = ^minl^ (2.42) 

where the minimum wavelength is 

^min ~ c/ fmai (2.43) 

and 

fmax = L'fo Hz. (2.44) 

Substituting equations (2.43) and (2.44) into equation (2.42) yields 

dr =dy = c/{2L'fo) (2.45) 

where L' is defined by equation (2.2) and fo is defined by equation (2.6). Substi- 
tuting equation (2.45) into equation (2.41) yields 



i = 0, il, i2, . . . 

d^g{q,m) ^ i:(qf L')mUo{q)7r + 2in, m = —M\ . . . .,0, . . . , M' (2.46) 

q = —L ', . . . , 0, . . . , L' . 
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By letting i = 0 in equation (2.46), the following rough estimate of direction cosine 
Uo at harmonic q and element m is obtained; 

Uo(q,m) = ± 9 7^0, m = ±1. (2.47) 

qmw 

For purposes of this thesis, the rough estimates of Uo{q) used were obtained by 
averaging the results obtained by evaluating equation (2.47) at elements m = ±1, 
that is, 

Uo(9) = 0.5(uo( 9, l) + Uo(?, -1)], q^O. (2.48) 



The most difficult case for the algorithm is to locate the highest harmonic, 
that is, q = ±X', when it exists at endfire (do = 90°) relative to the planar array. 
Let us assume that the value of the direction cosine at the highest harmonic is 
equal to one. Therefore, if ^ = +L\ equation (2.46) reduces to 



m) = mTT + 2i7r, 



i = 0,±1,±2,... 
m — —M ', . . . , 0, . . . , M' ’ 



(2.49) 



Evaluating this relationship at element m = ±1 to obtain the rough estimate of 
Uo{<l) yields 

C,(L',±1) = ±7r + 2i7r, i = 0, ±1, ±2, . . . , (2.50) 



and since the “wrapped” steady-state value must be in the interval [— tt, tt], i must 
be equal to zero. Equation (2.50) is reduced to 



C(i^',±l) = ±7r. 



(2.51) 



If the plus sign is chosen in equation (2.47), and m = ±1 and q = L' , then 
evaluation of equation (2.47) yields 

u,{L\±l) = eZ{L',±l)/±n). (2.52) 
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Upon substituting equation (2.51) into equation (2.52) we obtain 



Uo{L\ ± 1 ) = 1 



(2.53) 



which is the expected result in the absence of noise since the assumption was made 
that the value of the direction cosine was equeJ to one at the highest harmonic. 
Evaluating equation (2.49) at a different element (i.e., m ^ ±1) would yield i ^ 0 
and, as a result, an incorrect rough estimate of uq{L') = 1. This validates the 
use of elements m = ±1 to obtain the estimate of the direction cosines (equation 
(2.47)). 

The third step is to generate rough estimates of the “unwrapped” steady- 
state phase weights, B“^(q,m), by replacing Uo(q) in equation (2.46) with Uo(?) 
from equation (2.48), that is 



= ±(g/L’)mUo(q)7r, 



q — —L', . . . ,0, . . . , L' 
m = -A/',...,0,...,M'' 



(2.54) 



Therefore, if 



iir < 






<{i + 2)7t, 



t — 1 , 3 , 5 ,... 
q = —L', . . . , 0, . . . , 
m = —M', . . . ,0, . . . , M' 



(2.55) 



then 






i = 1,3,5, .. . 

5 = -L',...,0,...,L'(2.56) 
m = — M', . . . , 0, . . . , M' 
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where ] is the sign function. Note that the right-hand side of equation (2.56) 
is simply the noise corrupted version of the right-hand side of equation (2.35). 

The method of least-squares is then used to fit a straight line to the unwrapped 
phase weights 0^^{q, m) as a function of element number m for each harmonic q. The 
phase weights 0^g{q,m) were computed using equation (2.56). The least-squares 
slope at harmonic q is given by 



M' M' 

^Ls(q)= 

m — — M\ m = — 



( 2 . 57 ) 



Substituting equation (2.45) into equation (2.35) yields 



^“ 5(97 = ±{q/ L')Tmo{q)m, 



q = -L',...,0,...,L' 
m = — M', . . . , 0, . . . , M' 



(2.58) 



in the absence of noise. Therefore, since in the presence of noise 



^ ^Ls(9)»72 -f his{q), 



q — —L\ . . . , 0, . . . , L' 
m = — M', . . . , 0, . . . , M' 



(2.59) 



where 

1 

hs(q) = Tf E ^ss( 9 .">) (2.60) 

m = -M' 

is the least-squares “y intercept” at harmonic y, comparing equations (2.58) and 
(2.59) yields 

±«“(,) = (L7(,,r))lts(s). 9 5^0 ( 2 . 61 ) 



which is the least-squares estimate of direction cosine Uq at harmonic q. 

In a similar derivation, it can be shown that the least-squares estimate of the 
direction cosine Vo{q) at harmonic q is given by 



± 0 “ ( 5 ) =[i 7 (?ir)| iis(,), , ,^0 ( 2 . 62 ) 
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where now 



N' N’ 

SLs(q)= ^ ^ 

n=^N* n^—N* 

and, if 



(2.63) 



ITT < 



ks(q>^) 



< (i + 2 ) 7 t , 



i = 1,3,5, .. . 



(2.64) 



then 



i = 1 , 3 , 5 , . . . 

= Ks(q^^) + (* + l)7T, 9 = -L',...,0,...,I'(2.65) 



where 



= ±(9/i^')”^o(g)7r, 



9 — —L'l . . . , 0, . . . , L' 



n = -iV',...,0,...,7V' 

u„(g) = 0.5 [ 60 ( 9 , 1 ) + Uo(9,-l)], 9^0 



( 2 . 66 ) 

(2.67) 



and 

i^o(9,«) = ±-— <^«(9,«), 9^0, n = ±1. (2.68) 

^HTT 

Using the leeist-squares estimates of the direction cosines u^^(q) and Vo^{q)t 
that is, equations (2.61) and (2.62), respectively, yield estimates of the depression 
and bearing angles, ^ 0 ( 9 ) and V’o(9), respectively, using equations (2.33) and (2.34). 
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m. RESULTS 



Significant computer simulation results from four unique test cases are pre- 
sented in this chapter. The test cases were designed to test the algorithm’s multiple 
broadband target localization capability, its full angular coverage capability, and 
its angular resolution as a function of the input SNR at a single element in the 
array, sampling rate, harmonic number, and the number of iterations. 

The test cases were comprised of the following; 

(1) Case 1 - a single broadband target located at broadside relative to the planar 
array; 

(2) Case 2 - three broadband targets located at random positions; 

(3) Case 3 - a single broadband target located at endfire relative to the planar 
array; and 

(4) Case ,4 - three broadband targets, two of which share a common harmonic. 

The simulation results of these four test cases are based upon processing the 
output electrical signals from a 7 x 7 planar array of equally spaced hydrophones. 
The acoustic field incident upon the planar array was, in general, the sum of several 
plane-wave fields travelling in different directions. Each plane-wave field consisted 
of an arbitreiry number of hcirmonics (spectral lines) emanating from one of the 
broadband sound sources. As a result, the output electrical signal at each element 
of the planair Jirray was composed of an arbitrary number K of harmonics, all with 
identical amplitudes of unity. 

The fundaunental frequency for all test cases was chosen to be /„ = 1000 
Hz. Via equation (2.6), the fundamental period at each element of the array was 
1 millisecond. The sampling parameter, S, was set equal to 2, 4, and 6, which 
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corresponded to a sampling frequency equal to the Nyquist Rate, twice the Nyquist 
Rate, and three times the Nyquist Rate, respectively. The number of time samples 
taken per element of the array is given by 

L = SK + \ (3.1) 

where L is the total number of samples taken per element, K is the total number 
of harmonics present in the signal, and S is the sampling parameter. 

Baseline, or “no noise” test case results, were the first generated to ensure that 
the algorithm was working properly in a noise-free environment. Figure 3.1 depicts 
the noise free, time domain, received signed for case 2. For each baseline test case, 
and for a given value of the sampling parameter, bearing and depression angle 
estimation errors (measured in degrees) were obtained by running the computer 
simulation once and allowing the modified complex LMS algorithm 100 iterations. 
Case 1 baseline test case results showed zero degree estimation errors whereas case 
2 and case 3 baseline results showed estimation errors of less than 0.1 degrees. Case 
4 baseline results will be discussed later in this chapter. 

Following compilation of the baseline results, identical test cases were run 
using additive, wide-sense stationary, zero mean, white, guassian noise samples to 
corrupt the time samples of the received signal. Figure 3.2 depicts the time domain, 
received signal for case 2 for a SNR = 0 dB. For each test case, and for a given 
value of the sampling pairameter and input SNR at a single element of the array, 
average bearing and depression angle estimation errors were obtained by running 
the computer simulation 100 times per SNR value and allowing the complex LMS 
algorithm 100 iterations per run. In all test cases, allowing the adaptive algorithm 
more them 100 iterations did not decrease the estimation errors. 
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REAL RECEIVED SIGNAL AT ELEMENT (1,1) 
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Figure 3.1. Real Received Signal at Element (1,1), for Case 2 with No 
Noise Present. 
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REAL RECEIVED SIGNAL AT ELEMENT (1,1) 



O 




dB. 
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Case 1 placed a single target at broadside relative to the pleinar array ( i.e., 
9o = 0®). This was considered the simplest case for the algorithm since it is 
at broadside that the far-field beam pattern has its narrowest beamwidth and, 
as a result, the algorithm should provide good angular resolution [Ref. 8: pp. 
62-66]. The general plane- wave field radiated by the target consisted o{ K = 6 
harmonics. Therefore, with K = 6 and 5 = 2,4, and 6, only L = 13, 25, and 37 
time samples, respectively, were taken per element in the array (equation (3.1)). 
Figure 3.3 presents the average estimation error of the depression angle versus the 
sampling parameter 5 for the 0 dB SNR case. It can be seen that as the value of 
S increases for a given harmonic, the magnitude of the estimation error decreases. 
If S is held constant, it can be seen that as the harmonic number (q) increases, the 
magnitude of the estimation error decreases. These two trends were expected. As 
the value of S increases, more time samples are being processed and, as a result, 
the noise component of the received signal should tend to average out to be zero. 
Since more data is available for evaluation, a better estimate of the deterministic 
signal results. An increase in the harmonic number, q, represents an increase in the 
frequency and, as a result, the beamwidth of the far-field beam pattern decreases 
[Ref. 8; pp. 58-62]. The decreased beamwidth increases the cingular resolution of 
the algorithm and, as a result, decreases the estimation error. 

Figure 3.4 presents the results for case 1 at 9 dB SNR. The trends which 
were present in the 0 dB case are again apparent here. It also should be noted 
that as the SNR value increases (i.e., from 0 dB to 9 dB), the magnitude of the 
estimation error decreases. 

Results for bearing eingle estimation error in the broadside case are irrelevant 
since the target is directly above the array, and the bearing angle has no meaning. 
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Figure 3.3. Average Depression Angle Estimation Error vs Sampling 
Parameter 5, Case 1: SNR = 0 dB. I — 100. 
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Figure 3.4. Average Depression Angle Estimation Error vs Sampling 
Parameter Case 1: SNR = 0 dB. / = 100. 
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Ceise 2 placed three broadband targets at random locations. The general plane- 
wave field radiated by eaw:h of the targets contained two unique harmonics. Target 

1 was located at [Og = 49°, = 38°) and radiated harmonics 1 and 6. Target 

2 waa located at {9g = 5,xl>o = 137°) and radiated harmonics 2 and 5. Target 

3 was located at [Og = 77°, = 307°) and radiated harmonics 3 and 4. Since 

there are three incident plane-wave fields each containing two unique harmonics, 

the output electrical signal from each element in the array exhibits a total of = 6 
harmonics. As before, equation (3.1) dictates that only L = 13,25, and 37 time 
samples are to be taken per element of the array when S = 2,4, and 6, respectively. 
Figure 3.5 presents case 2 results for the depression angle estimation error with a 
0 dB SNR. In general, as S increases, the estimation error decreases for a given 
harmonic. When comparing the harmonics associated with a particular target, 
the highest harmonic usually is associated with the lessor of the estimation errors. 
These general results can be explained using the same arguments presented for case 
1 . 

Figure 3.6 presents the case 2 bearing angle estimation errors for 0 dB SNR. 
Figures 3.7 and 3.8 present the case 2 estimation errors for the depression and 
bearing angle, respectively, for 9 dB SNR. The same general trends already in- 
troduced are apparent in Figures 3.5 through 3.8. As in case 1, as the SNR value 
increased, the magnitude of the estimation errors decreased. 

Case 3 placed a single broadband target at endfire relative to the planar array 
{6g = 90°, 00 = 90°). This czise was considered the most difficult for the algorithm 
since it is at endfire that the feir-field beam pattern beamwidth is its broadest [Ref. 
8; pp. 62-66]. The general plane-wave field radiated by the target consisted of ^ = 5 
harmonics. Therefore, with 5 = 2,4, and 6, only L = 11,21, and 31 time samples 
were taken, respectively, at each element of the eirray (equation (3.1)). Figures 
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Figure 3.5. Average Depression Angle Estimation Error vs Sampling 
Parameter S, Case 2: SNR = 0 Db. / = 100. 
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nVERflGE ESTINfiTION ERROR FOR (OEGS) 
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Figure 3.6. Average Bearing Angle Estimation Error vs Sampling Pa^ 
rameter 5, Case 2: SNR = 0 dB I = 100. 
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Figure 3.7. Average Depression Angle Estimation Error vs Sampling 
Parameter Sy Case 2: SNR = 0 dB I = 100. 
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nVERFlGE ESTimilON ERROR FOR (DEGS) 
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Figure 3.8. Average Bearing Angle Estimation Error vs Sampling Pa- 
rameter 5, Case 2: SNR = 9 dB / = 100. 
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3.9 and 3.10 illustrate the average depression eind bearing angle estimation errors, 
respectively, for 0 dB SNR. Figures 3.11 and 3.12 illustrate the same information 
for 9 dB SNR. The estimation errors for the highest harmonic (q = 5) are not 
shown for any of the case 3 results since the magnitude of these errors was between 
70 and 90 degrees. This poor performance for the highest harmonic at endfire can 
be explained theoretically (see Chapter 2 - phase unwrapping). 

Case 4 corresponded to three broadband targets being present. The general 
plane-wave field radiated by each target was composed of two harmonics. In this 
case, two of the targets shared a common spectral line, namely, the fundamental 
frequency {q = 1). Target 1 was located at (6^ = 45®,t/>^ = 0°) and radiated 

harmonics 1 and 2. Target 2 was located at {9^ = 45°,V’„ = 180") and radiated 

harmonics 1 and 5. Target 3 was located at {9^ = 33°,V’q = 47°) and radiated 

harmonics 3 and 4. A total of AT = 5 harmonics were present in the output 

electrical signals, and, with 5 = 2, 4, and 6; only L = 11,21, and 31 time samples, 
respectively, were talcen at each element of the array (equation (3.1)). For this 
case, only no noise results were compiled, with the computer simulation running 
once and the complex LMS adaptive algorithm allowed 100 iterations. 

For harmonics q = 2,3,4, and 5 contained in the output electrical signal, 
the algorithm correctly identified the location of the targets with zero degrees 
estimation error. However, for the shared harmonic {q = 1), the algorithm located 
a false fourth target, which was located exactly between the actual locations of the 
targets whose radiated general plane-wave fields contained the shared harmonic, 
namely, targets 1 and 2. 

Tabular numerical data for test cases 1 through 4 is contained in The Ap- 
pendix. 
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Figure 3.9. Average Depression Angle Estimation Error vs Sampling 
Parameter 5, Case 3: SNR = 0 dB / = 100. 
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Figure 3.10. Average Bearing Angle Estimation Error vs Sampling Pa- 
rameter 5, Case 3: SNR — 0 dB I = 100. 
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Figure 3.11. Average Depression Angle Estimation Error vs Sampling 
Parameter 5, Case 3: SNR = 0 dB I — 100. 
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Figure 3.12. Average Bearing Angle Estimation Error vs Sampling Pa- 
rameter Sy Case 3: SNR — 9 dB I = 100. 
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test ccise 2 because it Weis felt that this ceise presented the most realistic real world 
situation. First, a minimum number of iterations for the algorithm to produce the 
same results as presented in this chapter (when 100 iterations of the algorithm was 
allowed) Wcis found. The minimum number was found to be 10 iterations. Using 
this information, the step size pajameter, n, (which had been a constant) was 
varied using a step function. For the first 10 iterations of the algorithm, n was held 
constant. This allowed the algorithm to reach its minimum value with a constant 
H. The value for n was then stepped down to one-tenth its original value for an 
additional 90 iterations. No decrease in estimation error resulted. The parameter 
was then stepped down to one one-hundredth, and then one one-thousandth, of 
its original value, with no change in the estimation error resulting. Additional 
attempts to model the step size parameter as a decreasing exponential yielded no 
change in the estimation error. 
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IV. CONCLUSIONS AND RECOMMENDATIONS 



The goal of this thesis was to evaluate a frequency domain adaptive beamform- 
ing algorithm which was developed by Ziomek and Chan [Ref. 5]. The multiple 
broadband target localization capability and the full angular coverage capability 
of the algorithm were of particular interest. 

Chapter 3 presented the results of four test cases which were designed to 
test the algorithm for several different capabilities. Severail conclusions concerning 
the performance of the frequency domain adaptive beamforming algorithm can be 
drawn from these results. Among these conclusions is the fact that the algorithm 
performs well in all test cases, yielding what would be considered acceptable esti- 
mation errors for the target localization problem. It should also be noted here that 
the number of time samples taken per element of the array never exceeded 37 time 
samples, an amount which represents a small number of data points. 

The algorithm does exhibit multiple broadband target localization capability. 
The test case results show that multiple targets can be localized if their radiated 
acoustic plane-wave fields contain at least one unique spectral line. The separation 
distance of the spectral lines can be controlled by selecting an appropriate value for 
the FFT bin spacing. The algorithm has also demonstrated its ability to localize 
targets both at endfire and broadside positions relative to the planar array, thus 
demonstrating the algorithm’s full cingular coverage capability. 

Several general trends were apparent in the results. The first of these trends is 
that as the sampling parameter S increased, the magnitude of the estimation error 
decreased. This wm due to a corresponding increase in the amount of time samples 
being processed. The second trend observed was that as the harmonic number q 
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increased, the magnitude of the estimation error, in general, decreased. This trend 
can be explained by the fact that as the harmonic number q increases (equivalent to 
an increcise in the frequency), the far-field beam pattern beamwidth decreases which 
results in better angular resolution. The third trend observed was that as the SNR 
increased (the output electrical signal becoming less noisy), the magnitude of the 
estimation error decreased. It was also found that the estimation errors presented 
in this thesis axe obtainable after ten iterations of the adaptive algorithm, and an 
increase in the number of iterations did not decrease the estimation error. 

In the course of this investigation, several possible areas for future research 
presented themselves: 

• further study of the effects of varying the step-size parameter, //, in an attempt 
to decrease the magnitude of the estimation error, 

• investigation of a noise reduction system prior to processing the output elec- 
trical signals, and 

• application of other spectral analysis techniques (i.e., autoregressive, maxi- 
mum entropy, maximum likelihood, etc.) to produce frequency spectra. 
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APPENDIX 



The Appendix presents the tabular numerical data from which the graphical 
results of Chapter 3 were produced. The numerical data for the graphical results 
of Figures 3.3 through 3.12 is presented in Tables A.l through A. 10, respectively. 
Tables A. 11 and A. 12 presents the numerical data for Case 4. There is no corre- 
sponding graphical result for Case 4. 



TABLE A.l. NUMERICAL DATA CORRESPONDING TO FIGURE 3.3. 



Case 1: One target located at broadside relative to planar 
array; six harmonics present in output electrical signal; 
Sampling parameter S = 2, 4, and 6; 1 00 iterations; SNR = OdB 



HARMONIC LOCATION 

q (DEG) (DEG) 

1 0 270 

2 0 270 

3 0 270 

4 0 270 

5 0 270 

6 0 270 



DEPRESSION ANGLE 6 ^ 
ESTIMATION ERROR (DEG) 



S=2 


5=4 


5=6 


-7,388 


-4.732 


-4.124 


-3.287 


-2.760 


-1.975 


-2.542 


-1.654 


-1.374 


-2.506 


-1.223 


-1.099 


-1.312 


-0.962 


-0.784 


-1.086 


-0.892 


-0.660 
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TABLE A.2. NUMERICAL DATA CORRESPONDING TO FIGURE 3.4. 



Case 1; One target located at broadside relative to planar 
array; six harmonics present in output electrical signal, 
Sampling parameter S = 2, 4, and 6; 1 00 iterations; SNR = 9dB 



HARMONIC LOCATION 

q (DEG) (DEG) 

1 0 270 

2 0 270 

3 0 270 

4 0 270 

5 0 270 

6 0 270 



DEPRESSION ANGLE 
ESTIMATION ERROR (DEG) 



S=2 


S=4 


S=6 


-2.312 


-1.796 


-1.456 


-1.191 


-0.940 


-0.686 


-0.769 


-0.582 


-0.449 


-0.594 


-0.446 


-0.365 


-0.479 


-0.326 


-0.267 


-0.410 


-0.316 


-0.214 



TABLE A.3. NUMERICAL DATA CORRESPONDING TO FIGURE 3.5. 



Case 2: Multiple targets located at random positions; 

six harmonics present in output electrical signal; 

Sampling parameter S = 2, 4, and 6; 100 iterations; SNR = OdB 



HARMONIC 


LOCATION 


DEPRESSION ANGLE 6 ^ 




e. 




ESTIMATION ERROR (DEG) 


q 


(DEG) 


(DEG) 


S=2 


S=4 


S=6 


1 


49 


38 


-2.789 


-0.708 


-0.285 


2 


5 


137 


-0.686 


-0.450 


-0.522 


3 


77 


307 


+ 7.138 


+ 3.601 


+ 1.075 


4 


77 


307 


+ 1.818 


+ 1.085 


-0.521 


5 


5 


137 


-0.599 


-0.164 


-0.010 


6 


49 


38 


-2.452 


-2.752 


-0.383 
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TABLE A.4. NUMERICAL DATA CORRESPONDING TO FIGURE 3.6. 



Case 2: Multiple targets located at random positions; 

six harmonics present in output electrical signal; 

Sampling parameter S = 2, 4, and 6; 1 00 iterations; SNR = OdB 



HRRMONIC 


LOCRTION 


BERRING RNGLE x\> 

’ O 




e<, 


4^, 


ESTIMATION ERROR (DEG) 


q 


(DEG) 


(DEG) 


S=2 


S=4 S=6 


1 


49 


38 


-0.261 


+ 0.236 +0.001 


2 


5 


137 


+ 0.270 


-0.043 -1.337 


3 


77 


307 


+ 43.268 


+ 27.662 +14.599 


4 


77 


307 


+ 9.620 


+ 3.789 +1.100 


5 


5 


137 


+ 0.168 


+ 1.122 +0.702 


6 


49 


38 


-0.420 


-3.220 +0.907 



TABLE A.5. NUMERICAL DATA CORRESPONDING TO FIGURE 3.7. 



Case 2: Multiple targets located at random positions; 

six harmonics present in output electrical signal; 

Sampling parameter S = 2, 4, and 6; 100 iterations; SNR = 9dB 



HRRMONIC 


LOCRTION 


DEPRESSION RNGLE 0, 




©0 


4^0 


ESTIMATION ERROR (DEG) 


q 


(DEG) 


(DEG) 


S=2 


S=4 


S=6 


1 


49 


38 


-0.338 


-0.1 13 


-0.080 


2 


5 


137 


-0.053 


-0.073 


-0.050 


3 


77 


307 


-0.425 


+ 0.280 


-0.192 


4 


77 


307 


+ 1.383 


+ 0.372 


+ 1.350 


5 


5 


137 


+ 0.017 


-0.013 


-0.003 


6 


49 


38 


-0.447 


+ 0.059 


+ 0.022 
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TABLE A.6. NUMERICAL DATA CORRESPONDING TO FIGURE 3.8. 



Case 2: Multiple targets located at random positions; 

six harmonics present in output electrical signal, 

Sampling parameter S = 2, 4, and 6; 1 00 iterations; SNR = 9dB 



HARMONIC 


LOCATION 


BEARING ANGLE 4^^ 




©o 




ESTIMATION ERROR (DEG) 


q 


(DEG) 


(DEG) 


S=2 


S=4 


S=6 


1 


49 


38 


■^0.235 


+0.048 


-0. 1 45 


2 


5 


137 


M.288 


+ 0.263 


+0.725 


3 


77 


307 


-0.346 


+0.037 


-0.036 


4 


77 


307 


-1.221 


-0.676 


-1.298 


5 


5 


137 


+0.648 


-0. 1 48 


-0.538 


6 


49 


38 


+0.058 


+0.1 12 


+0.105 



TABLE A.7. NUMERICAL DATA CORRESPONDING TO FIGURE 3.9. 



Case 3: One target located at endfire relative to planar array; 

five harmonics present in output electrical signal; 

Sampling parameter S = 2, 4, and 6; 1 00 iterations; SNR = OdB 



HARMONIC 


LOCATION 


DEPRESSION ANGLE 6, 






'i'o 


ESTIMATION ERROR (DEG) 


q 


(DEG) 


(DEG) 


S=2 


S=4 


5=6 


1 


90 


90 


+7.308 


+ 7.51 1 


+ 7.359 


2 


90 


90 


+ 16.978 


+ 15.058 


+ 8.046 


3 


90 


90 


+ 7.594 


+ 6.1 19 


+4.723 


4 


90 


90 


+ 13.749 


+7.554 


+2.926 


5 


90 


90 


+72.572 


+ 68.769 


+ 70.851 
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TABLE A.8. NUMERICAL DATA CORRESPONDING TO FIGURE 3.10. 



Case 3: One target located at endfire relative to planar array; 

five harmonics present in output electrical signal; 

Sampling parameter S = 2, 4, and 6; 1 00 iterations; SNR = OdB 



HARMONIC 


LOCATION 


BEARING ANGLE 


4^. 




0, 


4', 


ESTIMATION ERROR (DEG) 


q 


(DEG) 


(DEG) 


S=2 


S=4 


S=6 


1 


90 


90 


-0.547 


-0.373 + 


0.213 


2 


90 


90 


-19.202 


-18.163 - 


■7.252 


3 


90 


90 


+0.168 


-0.038 - 


•0.094 


4 


90 


90 


-3.171 


-0.715 - 


•0.041 


5 


90 


90 


-69.481 


-89.808 - 


•88.562 



TABLE A.9. NUMERICAL DATA CORRESPONDING TO FIGURE 3.11 



Case 3; One target located at endfire relative to planar array; 

five harmonics present in output electrical signal; 

Sampling parameter S = 2, 4, and 6; 1 00 iterations; SNR = 9dB 



HRRMONIC LOCATION DEPRESSION ANGLE 0^ 

0. 4^ ESTIMATION ERROR (DEG) 



q 


(DEG) 


(DEG) 


S=2 


S=4 


S=6 


1 


90 


90 


+5.301 


+ 4.994 


+ 4.677 


2 


90 


90 


+3.523 


+ 3.175 


+ 2.932 


3 


90 


90 


+ 3.744 


+ 2.827 


+ 2.701 


4 


90 


90 


+ 2.483 


+ 2.592 


+2.265 


5 


90 


90 


+68.938 


+ 66.092 


+ 70.858 
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TABLE A.IO. NUMERICAL DATA CORRESPONDING TO FIGURE 3.12. 



Case 3: One target located at endfire relative to planar array, 
five harmonics present in output electrical signal; 

Sampling parameter 5 = 2, 4, and 6; 1 00 iterations; SNR = 9dB 



HARMONIC 


LOCATION 


BEARING ANGLE 


4^0 




6o 


^0 


ESTIMATION ERROR (DEG) 


q 


(DEG) 


(DEG) 


S=2 


S=4 


S=6 


1 


90 


90 


+0.333 


+ 0.097 - 


0.105 


2 


90 


90 


-0.030 


-0.123 + 


0.088 


3 


90 


90 


-0.090 


-0.045 + 


0.019 


4 


90 


90 


+0.039 


+ 0,009 - 


0.021 


5 


90 


90 


-108.012 


-90.250 - 


93.038 



TABLE A.ll. NUMERICAL DATA FOR DEPRESSION ANGLE ESTIMATION 

ERRORS FOR CASE 4. 



Case 4: Multiple targets, targets one and tv/o share a common 
spectral line; five harmonics present in output electrical signal; 
Sampling parameter S = 2, 4, and 6; 1 00 iterations; no noise 



HARMONIC 


LOCATION 


DEPRESSION ANGLE 6, 




e. 




ESTIMATION ERROR (DEG) 


q 


(DEG) 


(DEG) 


S=2 


S=4 


S=6 


1 


45 


0 


+ 19.683 + 


19.683 


+ 19.683 


2 


45 


0 


0.0 


0.0 


0.0 


1 


45 


180 


+ 19.683 + 


19.683 


+ 19.683 


5 


45 


180 


0.0 


0.0 


0.0 


3 


33 


47 


0.0 


0.0 


0.0 


4 


33 


47 


0.0 


0.0 


0.0 
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TABLE A.12. NUMERICAL DATA FOR BEARING ANGLE ESTIMATION 

ERRORS FOR CASE 4. 



Case 4: Multiple targets; targets one and two share a common 
spectral line; five harmonics present in output electrical signal; 
Sampling parameter S = 2, 4, and 6; 1 00 iterations; no noise 



HRRMONIC 


LOCATION 


BEARING ANGLE 


4^. 




0. 


4'. 


ESTIMATION ERROR (DEG) 


q 


(DEG) 


(DEG) 


S=2 


S=4 


S=6 


1 


45 


0 


-180.00 


-180.00 - 


180.00 


2 


45 


0 


0.0 


0.0 


0.0 


1 


45 


180 


0.0 


0.0 


0.0 


5 


45 


180 


0.0 


0.0 


0.0 


3 


33 


47 


0.0 


0.0 


0.0 


4 


33 


47 


0.0 


0.0 


0.0 
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